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(57) ABSTRACT 

A mobile pill transmitter system that moves through or adja- 
cent to one or more organs in an animal’s body and that 
provides signals from which the pill’ s present location and/or 
present angular orientation can be determined. The system 
also provides signals from which the present roll angle of the 
pill, about a selected axis, can be determined. When the 
location coordinates and the roll angle of the pill are within 
selected ranges, an aperture on the pill container releases a 
selected chemical into or onto the body. Optionally, the pill as 
it moves also provides a sequence of visually perceptible 
images; the times for image formation may correspond to 
times at which the pill transmitter system location or image 
satisfies one or at least four different criteria. 

11 Claims, 3 Drawing Sheets 
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DETERMINATION AND APPLICATION OF 
LOCATION AND ANGULAR ORIENTATION 
OF A PILL TRANSMITTER WITHIN A BODY 

ORIGIN OF THE INVENTION 5 

This invention was made by an employee of the U.S. gov- 
ernment. The U.S. government has the right to make, use 
and/ or sell the invention described herein without payment of 
compensation, including but not limited to payment of roy- 10 
alties. 

FIELD OF THE INVENTION 

This invention relates to determination of location and 15 
angular orientation of a signal transmitting mobile device 
within a human or other animal body. 

BACKGROUND OF THE INVENTION 

20 

Medical imaging procedures now include, in some circum- 
stances, provision of a sequence of images, in video format, 
of the interior of an organ through which the imaging mecha- 
nism passes. The images provided are often life-like and 
contain considerable detail. However, as yet no such imaging 25 
mechanism has provided, at the same time, location coordi- 
nates and/or angular orientation coordinates, associated with 
this mechanism, indicating where the mechanism is presently 
located and/or presently oriented, when a particular image is 
made. This location and/or orientation information would be 30 
useful in associating a location and/or orientation with a 
particular image where a troublesome or pathological organ 
condition is indicated by the image. This location and/or 
orientation information would also be useful in determining, 
on a first pass or on a second pass, a location and/or orienta- 35 
tion where a particular medical treatment or pharmaceutical 
is to be delivered to the organ or a visibly perceptible image is 
to be made. 

What is needed is a method and system for determining 
location and angular orientation, exact and without iterations 40 
or approximations, of a mobile signal transmitting device, 
moving within a human or other animal body (“user”) at each 
of a sequence of times. Preferably, the system should con- 
sume minimal additional energy and should permit arbitrary 
separation of consecutive times in the time sequence(s). Pref- 45 
erably, the system should flexible to permit changes in the 
parameters to be monitored or examined. Preferably, the sys- 
tem should permit accurate determination of one or more 
times and corresponding locations for delivery of a drug or 
other chemical substance to the user. 50 

SUMMARY OF THE INVENTION 

These needs are met by the invention, which provides and 
applies an algorithm for exact determination of location coor- 55 
dinates and angular orientation coordinates for a mobile pill 
transmitter (“PT”) or other similar device that is introduced 
into and moves within a human or other animal body. A set of 
as many as eight nonlinear equations is developed and 
applied, relating propagation of a wireless signal between 60 
either of two, three or more or more transmitting antennas 
located on the PT to four or more non-coplanar receiving 
antennas located on a signal receiver appliance worn by the 
user. The equations are solved exactly, without approxima- 
tions or iterations, and are applied in several environments: 65 
(1) association of a visual image, transmitted by the PT at 
each of a second sequence of times, with a PT location and PT 


2 

angular orientation at that time; (2) determination of a posi- 
tion within the body at which a drug or chemical substance or 
other treatment is to be delivered to a selected portion of the 
body; (3) monitoring, after delivery, of the effect(s) of admin- 
istration of the treatment; and (4) determination of one or 
more positions within the body where provision and exami- 
nation of a finer scale image is warranted. 

The invention is applicable wherever a mobile pill trans- 
mitter is administered and followed within an organ of a 
human or other animal body. Provision of location and angu- 
lar orientation coordinates can be used wherever examination 
and/or therapeutic treatment and/or post -treatment examina- 
tion of the user is required. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 schematically illustrates relationships of a mobile 
pill transmitter and a signal receiving appliance operated 
according to embodiments of the invention. 

FIGS. 2A/2B are a flow chart of procedures for practicing 
the invention. 

DESCRIPTION OF BEST MODES OF THE 
INVENTION 

With reference to FIG. 1, a mobile pill transmitter {“PT”) 
11 is introduced into the body of a human or other animal 
(“user”), and the PT moves slowly within or adjacent to a 
body organ, for example, the stomach, small intestine, large 
intestine or other organ in the GI tract. First and second signal 
transmitters, 12-/ (i=l, 2), are provided at spaced apart loca- 
tions on the PT 11 with location coordinates (x0±a*Ax0, 
y0±a*Ay0,z0±a*Az0), where (x0,y0,z0) are the location coor- 
dinates of a specified center of the PT and a (0<a^ 1 ) is a 
selected real number (e.g., a=0.5 or a=l). Herein, it is 
assumed that a=l , although any other permitted choice can be 
used. The pill transmitter 11 is introduced into, and travels 
through a portion of, the user’s body to form and transmit 
visual images and/or other electronic signals, using a data 
antenna 13 that may coincide with, or be different from, one 
of the antennas 12-/ (i=l,2). 

A signal receiving appliance 14, spaced apart from the 
P[l], is provided with first, second, third and fourth spaced 
apart, signal receiving antennas, 15 -j (j=l, 2, 3, 4), which are 
non-coplanar so that the four antenna locations do not lie on 
a single plane. The first, second, third and fourth antennas 
have the respective location coordinates (xl,yl,zl), (x2,y2, 
z2), (x3,y3,z3) and (x4,y4,z4). The appliance 14 has one or 
more data-receiving antennas 16, which may be different 
from, or may coincide with, one or more of the four signal 
receiving antennas 15-y, to receive the data signals from the 
data antenna(s) 13. 

A location signal, transmitted from the transmitter antenna 
12-/ and received by the receiver antenna 15-j, will require a 
signal propagation time At(i,j) determined by the relation 

(xj-x0 ± AxO ) 2 + (yj-yO ±Ay0) 2 +(z/-z0±Az0) 2 = 

{c'A t(ij)} 2 , (1-i-j) 

where, for example, +Ax0 and -AxO correspond to i=l and 
i=2, respectively, c* is a representative velocity of signal 
propagation within the body, and At(i,j) is signal propagation 
time (measured or estimated) from transmitter no. i to 
receiver no. j. Squaring Eqs. (1-1-j) and (1-2-j) and forming 
the difference yields a new sequence of relations: 

4(xj-xO)AxO+4(yj-yO)AyO+4(zj-zO)AzO= 

c' 2 {At(lj) 2 -At(2j) 2 } (2-j) 
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Subtracting Eq. (2-1) from Eq. (2-2), Eq. (2-1) from Eq. (2-3) 
and Eq (2-1) from Eq. (2-4) provides three independent rela- 
tions, which are expressed in matrix form: 



AxO 


t Ar(l , 2) 2 - Ar(2, 2) 2 - Af(l, l) 2 + Ar(2, l) 2 } 

M 

AyO 

= (e- 2 /4) 

{Ar(l, 3) 2 - Ar(3, 3) 2 - Ar(l, l) 2 + Ar(3, l) 2 } 


A^O 


{Ar(l, 4) 2 - Ar(4, 4) 2 - Ar(l, l) 2 + Ar(4, l) 2 } 


AxO 


{Ar(l, l) 2 - Ar(2, l) 2 - Ar(l, 2) 2 + Ar(2, 2) 2 } 

AyO 

= (c’ 2 /4)M~ l 

{Ar(l, l) 2 - Ar(3, l) 2 - Ar(l, 3) 2 + Ar(3, 3) 2 } 

A^0 


t Ar( 1 , l) 2 - Ar(4, l) 2 - Ar(l, 4) 2 + Ar(4, 4) 2 } 


M = 


xl-xi 
x3 - xi 
x4-xi 


yl-yi 

y3-y\ 

yA-yl 


zX ~zi 
z3-zi 
z4 -zi 


(5) 15 


The precision matrix M has a non-zero determinant, because 20 
the four locations with coordinates are non-coplanar. Thus, 

Eq. (3) has a unique solution vector {AxO AyO, AzO}^. The 
square of the separation distance d between the two transmit- 
ting antenna locations is 25 

4{Ax0 2 +Ay0 2 +Az0 2 } =d 2 , (6) 


where d is a known, controllable physical distance. Equation 
(6) serves as a consistency relation to evaluate the accuracy of 
Eq. (4). 

One can also add Eq. (1-1 -i) to Eq. (1-2-j) and use Eq. (6) 
to obtain 


(xj - xO) 2 + (yj - yO ) 2 + (zj - zO) 2 + d 2 l4, 

= c’ 2 {Ar(l, jf + Af(2, jf}/2(j = 1, 2, 3, 4) 

and observe that 

xj-xO=(xj-xi)+(xi-xO)(j=i, 2, 3, 4). 

Using Eqs. (8-j), Eqs. (7-j) (j=2 , 3, 4) can be re-expressed as 

45 

(xi - xO) 2 + 2(xi - xO )(xj - xi) + (xj -xi) 2 + (9- j) 

(yi - yO) 2 + 2(yi - yO)(yj - yi) + (yj - yif + 

(zi - zO) 2 + 2(zi - zO)(zl - zi) + (zj - zi) 2 + d 2 / 4 = 50 

c ,2 {Ar(l, j) 2 + Ar(2, j) 2 }l2. 


(7-j) 


35 


40 


(8-j) 


Subtracting Eq. (9-1) from Eq. (9-j), for each of j=2, 3, 4, 55 
yields three linear relations in the remaining unknowns (j=2, 

3, 4): 


-continued 


xl -x£) 


{Ar(l, 2) 2 - Ar(2, 2) 2 - At(i, l) 2 + At(2, i) 2 } 

yl-yO 

= c ,2 M~ l 

{Ar(l, 3) 2 - At(3, 3) 2 - At(i, l) 2 + At(3, i) 2 } 

zi - zP 


{Ar(l, 4) 2 - At(4, 4) 2 - Af(l, l) 2 + Ar(4, l) 2 } 


The location coordinates (x0,y0,z0) of the center are then 
obtained by subtracting the known values of the antenna 
coordinates (x 1 ,y 1 , z 1 ) in Eqs .(12). 

Note that inversion of the precision matrix M is involved in 
each of the solution component triples {AxO AyO, Az0} rr and 
{x0,y0,z0}U The quantity S={det(M)} -1 is a measure of the 
sensitivity of these solutions to small errors in measurement 
and/or timing, which may arise from signal -to-noise prob- 
lems, from signal propagation problems, from location mea- 
surement errors and/or from timing errors. The sensitivity 
varies approximately monotonically with S, with S=0 being 
most desirable, so that, as the four receiver antenna locations 
approach a common plane (coplanar configuration), the sen- 
sitivity approaches an infinite value. For this reason, it is 
preferable to choose the location coordinates, (xl ,y 1 ,zl ), (x2, 
y2,z2), (x3,y3,z3) and (x4,y4,z4), for the four receiver anten- 
nas to maximize the value det(M)=l/S, consistent with the 
constraints imposed on the coordinates. 

The solutions of Eqs. (3) and (11) provide a determination 
of the location coordinates (x0,y0,z0) for the center of the PT 
device 11 and of the angular orientation (AxO, AyO, AzO) of this 
device relative to axes of a Cartesian coordinate system, 
which may include pitch and yaw angles of the device. 

Another quantity, roll angle 0 of the device about its lon- 
gitudinal axis AA, illustrated in FIG. 1, may be needed in 
some applications. To determine 0, it is convenient to provide 
a third signal transmitting antenna 12-3 on the PT, at a loca- 
tion that is not collinear with the locations of the other two 
signal transmitting antennas, 12-1 and 12 - 2 . 

If, as is likely, the PT is cylindrical and the roll angle 0 is to 
be estimated, it is preferable to locate the two transmitting 
antennas, 12-1 and 12 - 2 , at antipodal locations on the cylin- 
der, as illustrated in FIG. 1. With this configuration, the line 
L(l,2) connecting the two transmitting antennas, 12-1 and 
12 - 2 , is now oriented at an angle 

<t>=tan~ 1 {r0/ {Ax0 2 +Ay0 2 +Az0 2 } 1/2 } , (13) 

relative to the axis AA, Where rO is the radius of the 
cylinder. The third transmitting antenna 12-3 is assumed to be 
located at a known, selected rotational angle 0=03 (e.g., 
15°^03^9O°) relative to a longitudinal plane II (defined by 
roll angle 0=0) that contains the transmitting antennas 12-1 
and 12 - 2 . 

Where the cylinder or PT 11 has rolled or rotated by an 
angle 0 relative to its reference position (0=0), the distance 
Ar(i,j) from the first and second transmitting antennas 12-/ 
(i=l, 2) to the receiving antenna 15 -j is related to the elapsed 
time At(i,j) for the signal propagation by 


(*1 - *0 )(xj - *1) + (yi - y0)(yj - yl) + (zi - z0)(zj - zi) = (10- j) 

c ,2 {Ar(l, 2) 2 + Ar(2, 2) 2 } - Ar(l, l) 2 - Ar(2, l) 2 } 



xi — xO 


t Ar( 1 , 2) 2 - At(2, 2 ) 2 - At(i, i) 2 + At(2, l) 2 } 

M 

yl-yO 

= C ’ 2 

t Ar( 1 , 3) 2 - Af(3, 3) 2 - Ar(l, l) 2 + Ar(3, l) 2 } 


zi -z0 


t Ar( 1 , 4) 2 - At(4, 4) 2 - Ar(i, i) 2 + At(4, l) 2 } 


60 A r(i, j) 2 = {xj -x0-r0- cos (6 - 63)} 2 + (13- i- j) 

{yj -y0-r0- sin (9 - 93)} 2 + {zj - z0 + AzO) 2 , = c’ 2 At(i, j) 2 . 


65 and the distance from the third transmitting antenna 12-3 to 
the receiving antenna 35-/ is related to the elapsed time At(3,j) 
for the signal propagation by 
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Ar(3, j) 2 = {xj -xO-rO- cos(0)} 2 + (13-3- j) 

{yj - ^0 - rO ■ sin(0)} 2 + {zj - zO} 2 , = r 2 Ar(3, j) 2 . 


where the terms -i-AzO and -AzO in Eq. (13-j) correspond to 
i=l and i=2, respectively. 

The difference in distance squared Ar 2 between Eq. (13-i-j) 
and (13-3-j) becomes 

A r{i, 3 j/') 2 -Ar(3 y) 2 =±2(zz'-z0)Az0+Az0 2 -2r0(xz'-.x0) 

{cos(0-03)-cos 0}-2rO(yz'-_yO){siii(0-03)-sin 
Q}=c*{At(ij?-At(3j) 2 }. (14) 

Using multiple angle trigonometric identities, one verifies the 
relations 


cos{9 - 93) - cos 9 = (15) 

cos 0{cos03 - 1 } + sin# ■ sin 93 = (2(1 - cos 03)} 1/2 cos{0 - if/3), 

sin(0 - 93) - sin# = (16) 

sin 0{ cos 93 - 1} + cos 9- sin 93 = {2(1 - cos 03)} 1/2 sin{0- if/3), 

ift3 = tan -1 {sin 93 j (cos 93 - 1)} = {93 -n)l 2 (17) 


Equation (14) can be rewritten as 

Ar(z, j) 2 - Ar(3, j) 2 = (18) 

±2 (zi - zO)Az + A z 2 - 2r0{(xi - xO) 2 + { yi - }>0) 2 } 1/2 
{2(1 - cos 03)} 1/2 cos{0 - if/3 - if/) = c ,2 {Ar(z, j) 2 - At{3, j) 2 }, 

if/ = tan -1 {(_yz - y0)/{xi - *0)}. (19) 


All quantities in Eq. (18) are known except the roll angle 0, 
which will vary as the PT 1 1 moves along a path in or adjacent 
to an organ in the body. The quantity 0-if 3-ip in Eq. (1 8) can 
be estimated to within a multiplicative factor (±1) by mea- 
surement of the propagation time difference Ar(i,j) 2 -Ar(3,j) 2 

A dispensing aperture for a chemical (pharmaceutical or 
other substance), can be initially located at the initial rota- 
tional angle 0=0, can be allowed to move rotationally with the 
roll angle 0, and can be caused to dispense the chemical when 
the aperture is located adjacent to, and is directed at, a target 
section of an organ that is to be treated, at a target roll angle 
value. Given a knowledge of the roll angle 0 and of the 
location (x0,y0,z0) of the center of the PT 11, the chemical 
dispensation point is easily estimated. 

Optionally, the PT 11 also forms and transmits signals 
representing visually perceptible images at a sequence of 
times as the pill moves through the body, using the signal 
receiving antenna 16. In a first option, the antenna 16 is part of 
a signal transceiver that (1 A) determines when location coor- 
dinates of the pill lie within selected ranges, such as 

( ^xO(t)^xJ2,yfl ^y0{t)^y/2,zfl ^z0 {t)^zfl), (20) 

and (IB) when the pill location coordinates lie within these 
ranges, a signal is transmitted at the transceiver that causes the 
pill transmitter system to form and transmit a visually per- 
ceptible image. In a second option, the PT forms and trans- 
mits an image at two consecutive image transmission times, t l 
and t 2 , when 

{(x0(r 2 )-x0(r 1 )) 2 +0;0(r 2 )->;0(z i 1 )) 2 +(z0(r 2 )- 

zO(?i)) 2 } 1/2 =Ar(thr), (21) 


6 

where Ar is a selected threshold distance. In a third option, the 
PT forms and transmits an image at each of a selected 
sequence of times {t„}„, where 

l ? «+i-0=A*(thr), (22) 

and At(thr) is a selected separation time. In a fourth option, a 
numerical measure of image difference ID(t x ;t 2 ) is provided 
between an image formed at each of two times, t l and t 2 (>tj, 
and the image I(t 2 ) is transmitted when 

10 

IID(rl;r2)I^A/(thr), (23) 

where Al(thr) is a selected image difference threshold. Each 
of these image transmission options is made possible by 
j 5 availability of the location coordinates (x0(t),y0(t),z0(t)), or a 
related quantity. 

FIG. 2 is a flow chart covering procedures for practicing the 
invention. In step 21, the system receives, at the receiving 
antennas 15-y (j=l, 2, 3, 4 in FIG. 1), signals transmitted by 
20 the transmitting antennas 12-/ (i=l, 2) and measures the sig- 
nal propagation times At(i,j) from transmitter no. i to receiver 
no. j. Each transmitter may transmit a signal having a differ- 
ent frequency, indicium or other signal signature, or the two 
sets of signals (i=l, 2) may be transmitted in different, non- 
25 overlapping times intervals, with temporal separations 
At(sep)=0.001-5 sec, or more if desired. In steps 22 and 23, 
the system respectively computes or otherwise provides the 
quantities (c ,2 /4) { At(l ,j) 2 — At(2,j) 2 — At(l , 1) 2 +At(2, l) 2 } (j=2, 
3, 4) and the entries of the matrix M -1 , where M is the 
30 precision matrix, set forth in Eq. (5). In step 24, the system 
computes the angular orientation coordinates (AxO, Ay0,Az0) 
as set forth in Eq. (4), where c' is a measured or estimated 
velocity of signal propagation within and adjacent to the 
body. 

35 In step 25, the PT system computes or otherwise provides 
the quantities c ,2 { At(l ,2) 2 +At(2,2) 2 — At(l , 1 ) 2 — At(2, 1 ) 2 } - In 
step 26, the PT system computes the center location coordi- 
nates (x0,y0,z0) of the PT relative to a location (xl,yl,zl) of 
one of the receiving antennas, as set forth in Eq. (12). From 
40 these values, one can determine the center location and angu- 
lar orientation of the PT, relative to a Cartesian or other 
coordinate system and relative to the location (xl,yl,zl) of 
one of the receiving antennas. In step 27 (optional), the sys- 
tem computes the square of the transmitting antenna separa- 
45 tion distance, 4{Ax0 2 +Ay0 2 +Az0 2 }, and compares this with 
the square of the known separation distance d 2 , to evaluate the 
accuracy of the determinations in one or more of the steps 
21-25 

In step 28 (optional), c ,2 {At(i,j) 2 -At(3,j) 2 }-{Az0 2 ±2(zi- 
50 z0)Az0}, 2r0{(xi-x0) 2 +(yi-y0) 2 } 1/2 {2(l-cos 03)} 1/2 , ijA 
and if are computed or otherwise provided for at least one of 
i=l or 2. In step 29 (optional), the relation 


55 2r0{(j« - jcO) 2 + (yi - .y0) 2 } 1/2 {2(l - cos 03)} 1/2 cos{0 -if/3 -ip)} = d 1 8 ) 

-c ,2 {Ar(z, j) 2 - Ar(3, y) 2 } + {A^O 2 ± 2 {zi - z0)Az0} 


60 from Eq. (18) is used to estimate the quantities cos {0-if 3 -if ) 
and the roll angle 0, where 03, if 3 and if are quantities 
discussed in the preceding. In step 30 (optional), the system 
determines if the center location coordinates (x0,y0,z0) lie 
within selected ranges (xrl ^x0^xr2.yrl ^y()=yr2. 

65 zrl ^zO^zrs) and/or the roll angle 0 lies in a selected range 
(0rl^0^0r2) so that a selected chemical can be released 
through the aperture into or onto the body. If the answer to the 
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query in step 30 is “yes,” the PT system releases the selected 
chemical, in step 31 (optional). If the answer to the query in 
step 30 is “no,”, the system moves to the next step, if any. 

In step 32 (optional), the PT system determines if one (or 
more than one) of the selected criteria set forth in Eqs. (20), 5 
(21), (22) and (23) is satisfied. If the answer to the query in 
step 32 is “yes,” the PT system forms and transmits a visually 
perceptible image, in step 32 (optional). If the answer to the 
query in step 32 is “no,” the system takes no further action at 
that time, in step 34 (optional). 10 

The invention thus provides: (1) location coordinates and/ 
or angular orientation coordinates for the PT, as a function of 
time; (2) an estimation of roll angle for the PT, as a function 
of time; (3) an estimation of one or more times a selected 
chemical is appropriately delivered to an organ; and (4) deter- 1 5 
mination of a sequence of times for formation and transmis- 
sion of a signal representing an image sensed by the PT. 

What is claimed is: 

1. A method for determination of at least one location 20 
coordinate and at least one angular orientation coordinate of 
a mobile device that moves within or adjacent to at least one 
organ in a human or other animal body, the method compris- 
ing: 

(1 ) providing first and second spaced apart signal transmit- 25 
ting antennas, having index values i=l and i=2, respec- 
tively, on the device, spaced apart along a selected axis 
of the device, and transmitting a first sequence and a 
second sequence of distinguishable signals at the first 
and second transmitting antennas, respectively, where 30 
the first and second transmitting antennas have location 
coordinates (x0+Ax0,y0+Ay0,z0+Az0) and (xO-AxO, 
y0-Ay0,z0-Az0), respectively, a central location of the 
mobile device has location coordinates (x0,y0,z0), and 
the coordinate values xO, yO, zO, AxO, AyO, AzO are to be 35 
determined, and (x,y,z) are location coordinates in a 
selected Cartesian coordinate system; 

(2) receiving the first and second sequences of transmitted 
signals at least first, second, third and fourth spaced 
apart, non-coplanar signal receiving antennas, having 40 
index values j=l, 2, 3, 4, respectively, having substan- 
tially fixed and known location coordinates (xj,yj,zj), 
and being located adjacent to the body; 

(3) for each transmitted signal in each of the first and 
second sequences, providing an elapsed time of signal 45 
propagation At(i,j) between the first transmitting 
antenna (i=l) or second transmitting antenna (i=2) and 
the first, second, third or fourth signal receiving antenna 
(j=l, 2, 3, 4, respectively), where At(i,j) is expressed by 

a relation 


-continued 


' xl - *0 


' {Af(l, 2) 2 - Af(2, 2) 2 - Ar(l, l) 2 + Ar(2, l) 2 } ' 

o 

1 

= (c ,2 /4)M -1 

{Af(l, 3) 2 - Af(3, 3) 2 - Ar(l, l) 2 + Ar(3, l) 2 } 

. zl - zO 


{Af(l, 4) 2 - Af(4, 4) 2 - Ar(l, l) 2 + Ar(4, l) 2 } 
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(5) interpreting the quantities (xl-xO, yl-yO, zl-zO) as 
location coordinates of the mobile device central loca- 
tion relative to location coordinates (xl,yl,zl) of the 
first receiving antenna; and 

(6) when D={(AxO) 2 +(AyO) 2 +(AzO) 2 } 172 , is non-zero, 
interpreting the quantities (AxO/D, Ay 0/D, AzO/D), 
where as direction cosines with respect to x-, y- and 
z-Cartesian coordinate axes, corresponding to angular 
orientation of a line segment along the selected axis. 

2. The method of claim 1, further comprising evaluating 
accuracy of said angular orientation components by compar- 
ing a sum of squares of said angular orientation components, 
4{Ax0 2 +Ay0 2 +Az0 2 }, with a square d 2 of a known distance 
between two selected locations on said mobile device. 

3. The method of claim 1, further comprising estimating a 
roll angle 0 of said mobile device about a selected axis by a 
process comprising: 

selecting a location on said mobile device corresponding to 
a roll angle value, 0=03, about a selected axis for said 
mobile device; 

providing a third signal transmitting antenna, numbered 
i=3, on said mobile device at a location selected with 
reference to the roll angle value 0=03; 

providing a square At(3,j) 2 of elapsed time of signal propa- 
gation between the third signal transmitting antenna and 
at least one of said signal receiving antenna, numbered j , 

providing a first selected angle value xp =(0 3 — jc)/ 2 and a 
second selected angle value 

providing a value AzO for said angular orientation compo- 
nent and values, xi-xO, yi-yO and zi-zO, relating said 
location coordinate (x0,y0,z0) to corresponding loca- 
tion coordinate (xi,yi,zi) (1=1 or 2) of said signal receiv- 
ing antenna; and 

estimating a roll angle 0 of said mobile device at least one 
time from a relation 


2r0{{xi - xO) 2 + (yi - y0) 2 } 1/2 {2(! - cos 03)} 1/2 cos{03 - ij/3 - ijj) = 


c a At(i=l J) 2= (xj-x0+Ax0) 2 +(yj-y0+Ay0) 2 +(zj-z0+ 
AzO) 2 , 


—c ,2 {At{i, j f - Ar(3, y) 2 } + {AzO 2 ± 2 (zi - z0)Az0}, 


c’ 2 At(i=2j) 2 =(xj-x0-Ax0) 2 +(yj-y0-Ay0) 2 +(zj-z0- 

AzO) 2 , 55 

where c' is a representative value of propagation of a trans- 
mitted electromagnetic signal within the body; 

(4) estimating the values AxO, AyO, AzO, xO, yO and zO by 
the relations 60 


' A^0 ' 


' {Ar(l, l) 2 - Ar(2, l) 2 - Ar(l, 2) 2 + Ar(2, 2) 2 } ' 

AyO 

= (c ,2 /4)M -1 

{Ar(l, l) 2 - Ar(3, l) 2 - Ar(l, 3) 2 + Ar(3, 3) 2 } 

.AzO. 


{Af(l, l) 2 - Ar(4, l) 2 - Ar(l, 4) 2 + Ar(4, 4) 2 } 


where r0 is a radius of said mobile device, measured in a 
direction substantially perpendicular to the selected axis for 
said mobile device. 

4. The method of claim 3, further comprising choosing said 
angle xp to be x|j=tan _1 {(yi-y0)/(xi-x0)}. 

5. The method of claim 3, further comprising: 

providing said mobile device with a controllable aperture 

and with a source of a selected chemical associated with 
the aperture; 

estimating a time when (1) said location coordinates (xO, 
y0,z0) for said mobile device lie in selected ranges 
(xrl ^x0^xr2,yrl =y()^yr2,zrl ^z0^zr2) and (2) said 
roll angle 0 lies within a selected range, 0rl ^ 0^r2; and 
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when the conditions (1) and (2) are satisfied, causing said 
mobile device to release the selected chemical through 
the aperture. 

6. The method of claim 1, further comprising: 

causing said mobile device to form and transmit signals 
representing visually perceptible images at a sequence 
of two or more image formation times, when a selected 
criterion is satisfied. 

7. The method of claim 6, further comprising choosing said 
selected criterion to comprise: each image formation time 
corresponds to a time, t=t^ at which said location coordinates 
(XO(t^) ,y 0(ty) ,zO(ty)) lie within one or more selected ranges 
(xfl ^x0^xf2,yfl ^y0^^yf2,zfl ^z0^zf2). 

8. The method of claim 6, further comprising choosing said 
criterion to comprise: two consecutive image formation 
times, t=t^ and t=ty 2 , correspond to locations (xO^jyO^jZO 
(t^)) of said mobile device that satisfy 

{(x0(^)-^(/,) 2 + (><)(^)-^0(^) 2 + (z0(^)- 
z0 ({/ 2 ) 2 } 1/2 = Ar(thr), 

where Ar(thr) is a selected threshold distance. 


10 

9. The method of claim 6, further comprising choosing said 
criterion to comprise: two consecutive image formation time, 
tr=t^ and t=ij 2 satisfy 

l{/2 - {/ 1 !=A?(thr), 

where At(thr) is a selected threshold separation time. 

10. The method of claim 6, further comprising choosing 
said criterion to comprise: 

providing a numerical measure ID(t x ; t 2 ) of differences in 
; an image formed by said mobile device at times t=t 1 and 

t=t 2 (>ti); and 

the measure ID at two consecutive image formation times, 
andt^t^ satisfies 
ID(^)^ID(thr), 

15 where AID(thr) is a selected image difference threshold 
value. 

11. The method of claim 1, further comprising choosing at 
least one of said sets of coordinate values, (xl ,y l,zl), (x2,y2, 
z2), (x3,y3,z3) and (x4,y4,z4), to produce a magnitude of an 

20 inverse of a determinant, l/ldet(M)l , of said matrix M having 
a small positive value. 



